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Abstract 

An intense research on financial market microstructure is presently in progress. Continuous time 
random walks (CTRWs) are general models capable to capture the small-scale properties that high 
frequency data series show. The use of CTRW models in the analysis of financial problems is quite 
recent and their potentials have not been fully developed. Here we present two (closely related) 
applications of great interest in risk control. In the first place, we will review the problem of 
modelling the behaviour of the mean exit time (MET) of a process out of a given region of fixed 
size. The surveyed stochastic processes are the cumulative returns of asset prices. The link between 
the value of the MET and the timescale of the market fluctuations of a certain degree is crystal 
clear. In this sense, MET value may help, for instance, in deciding the optimal time horizon for 
the investment. The MET is, however, one among the statistics of a distribution of bigger interest: 
the survival probability (SP), the likelihood that after some lapse of time a process remains inside 
the given region without having crossed its boundaries. The final part of the article is devoted to 
the study of this quantity. Note that the use of SPs may outperform the standard "Value at Risk" 
(VaR) method for two reasons: we can consider other market dynamics than the limited Wiener 
process and, even in this case, a risk level derived from the SP will ensure (within the desired 
quintile) that the quoted value of the portfolio will not leave the safety zone. We present some 
preliminary theoretical and applied results concerning this topic. 
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I. INTRODUCTION 



The continuous time random walk (CTRW) formalism was introduced four decades ago 

n 

by Montroll and Weiss jlj, as a natural extension of ordinary random walks (RWs). In a (one 
dimensional) RW you can randomly move through a fixed grid either up or down, at regular 
time steps, whereas in a CTRW the size of the movements and specially the time lag between 
them are random. CTRWs have been successfully applied to a wide and diverse variety of 
physical phenomena over the years |2j: transport in random media, random networks, self- 
organized criticality, earthquake modelling; and recently also to finance 
llfli II, Q|. In this latter context, the efforts have been mostly focused on the statistical 



properties of the waiting time between successive transactions and the asset return at each 
transaction. Different studies in different markets are conceiving the idea that the empirical 
distributions of both random variables are compatible with an asymptotic fat tail behaviour 

Within the CTRW formalism we have recently investigated the mean exit time (MET) of 
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. In these articles we show 



asset prices out of a given region for financial time series 
that the MET follows a quadratic growth in terms of the interval width, both in small and 
large scales. We checked the persistence of this behaviour in time series from several markets, 
such the foreign exchange market, or the New York Stock Exchange (NYSE). The theoretical 
model used in these works was based on two-state chain Markovian processes. This model 
is able to both describe the quadratic scaling property observed for the MET and provide a 
mechanism that can incorporate asset peculiarities through return autocorrelations. 

One of the possible applications of the analysis of the MET in finance is in the field of 
risk control. There is a direct link between the value of the mean exit time out of a region, 
and the timescale of market fluctuations of a certain size. Therefore, its value may help, 
for instance, in deciding the minimal time horizon for an investment, the rotation rate of a 
portfolio, or even the value of stop-loss and stop-limit levels for a position. 

However, the mean exit time is only a statistic of a distribution with even bigger interest: 
the survival probability (SP), the probability that after some elapsed time a process remains 
inside the given region without having crossed its boundaries. This quantity may outperform 
the standard "Value at Risk" (VaR) method for two reasons: it could be based on market 
statistics different than the (unrealistic) Gaussian distribution, and it will ensure (within 
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the desired quintile) that the market value of the portfolio will not leave the safe zone. 

The paper is organized as follows. In Sect. |H]we discuss the MET within the CTRW 
formalism, under a meaningful set of simplifying assumptions. In Sect. IHII we relax some of 
the previous constrains in order to introduce some memory into the process. Section llVl is 
devoted to the SP, its properties and its connections with the MET. In Sect. El we show in 
a practical situation how SP can be used in risk control. Conclusions are drawn in Sect. IVI1 



II. EXTREME EVENTS WITHIN CTRW 



In the most common version of the CTRW formalism a given random process X(t) 
shows a series of random increments or jumps at random times • ■ • , £_i, to, ti, t2, ■ ■ ■ , t n , ■ ■ ■ 
remaining constant between these jumps. Therefore, after a given time interval r n = t n —t n -\, 
the process experiences a random increment AX n (r n ) = X(t n ) — X(t n _i) and the resulting 
trajectory consists of a series of steps as shown in Fig.^ Waiting times r n and random jumps 



AX„(r„) are described by their probability density functions (pdfs 



which we wi 



1 denote by 



l(j for a 



ip(r) and h(x) respectively. We refer the reader to Refs. 
more complete account of the CTRW formalism. 

In the present work we will show two applications of CTRWs to the study of extreme 
problems in financial time series. We will take as underlying random process X(t) the 
logarithmic price X(t n ) = ln(S(i n )), where S(t) is the stock price at time t. We first 
consider the problem of obtaining the mean exit time of X(t) out of a given interval [a, b], 
of width L. We assume that at certain reference time t , right after an event, the price has 
a known value X(to) = xq, xq G [a, b]. Let us focus our attention on a particular realization 
of the process and suppose that at certain time t n > to the process first leaves the interval 
— see Fig. ^ We call the lapse t n — to, the exit time out of the region [a, b] and we will 
denote it by t\ a u(xo)- This quantity is a random variable since it depends on the particular 
trajectory of X(t) chosen and the MET is simply the average T^(x ) = E[t[ afe ](a;o)]. 

The standard approach to exit time problems is based on the knowledge of the survival 
probability — see Sect. IIVI In general, this is a quite involved path ^]. However, if we 
assume that r n and AI„(r„) are independent and identically distributed (i.i.d.) random 
variables, described by a joint pdf p(x, r), 

p(x, r)dxdr = Probja; < AX n < x + dx; r < r n < r + dr}, 
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FIG. 1: A sample trajectory of the X{t) process along with the corresponding value of the random 
variable t[ a ^ (xq). 



it can be shown 11] that one can obtain the MET directly, without making use of the survival 
probability. In this framework the MET T^ a ^{xo) obeys the following integral equation: 



T[ a ,b]{xo) = E[t] + / h(x - x )T [aM (x)dx, 

J a 



(1) 



where E[r] is the mean waiting time between jumps. It is worth noticing that Eq. (P) is still 
valid even when r n and AX n are cross-correlated. In fact, in the case of an i.i.d. process 
the MET only depends on p(x, r) through its marginal pdfs ip(r) and h(x). 

We can illustrate the problem with a choice for h(x), based on the small-scale properties 



of the system, which results in the observed 
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quadratic growth in the MET. Let us 



introduce the following symmetrical two-state discrete model 

h(x) = - [5(x — c) + 5(x + c)] . 



cgr< 



(2) 



where c is the basic jump size. This choice for h(x) implies that the flat levels in every 
particular trajectory will be in a regular grid of size c centred at the starting point xq. It is 



worth noticing that this approach is also used in the context of option pricing, when the fair 
price of a derivative product is obtained by making use of the binomial trees methodology, 
where it is assumed that the stock price makes a jump up or down with some probability [14]. 
The solution of this problem, if we start from the middle of the interval, reads: 



T M (q + L/2) 
E[r] 



1 



L 



(3) 



If we consider a symmetric exponential function for the jump distribution instead: 

h(x) = Je" 7|a;| , (4) 



a very similar result is obtained ^j] : 

T M (q + L/2) = I 
E[t] 2 



1 + 



7 L 



(5) 



III. MEAN EXIT TIME FOR MARKOV-CHAIN MODELS 

In order to embrace also CTRWs with memory, we derived in [12] an integral equation for 
the MET when the jumps are Markovian. In particular, we focused on the case in which it 
is possible to neglect the influence of the past waiting time by assuming that the magnitude 
of the previous change carries all the relevant information. The equation in this case is: 



T [aM {x \AX ) = E [t\AX ] + f h{x- x \AX )T [aM (x\AX)dx, 

J a 



(6) 



with AX = x— Xq. Now the MET depends only on the marginal pdf of the return increments, 
h(x\AX ), and on the conditional expectation of the waiting time, E [t|AX ], which has to 
be evaluated through the marginal pdf, ^(t|AXq). In order to solve Eq. © and obtain 
explicit expressions for the MET we will use again a discrete two-state model: 

h{x\y) = 2c 5{x - c) + 2c 5{x + c), 

where r is the correlation between the magnitude of two consecutive jumps. The MET 
starting from the middle of the interval reads now: 

T [aM {a + L/2) 2r A L\ 1-r 



1 + r V 2c 



L 



E[r] 1 + r \ '2cy'l + r 

and, for large values of L/c, we recover the quadratic behaviour in the leading term: 

T M (a + L/2) 1-r 



EH 



1 / V 1 + 2c 
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IV. THE SURVIVAL PROBABILITY 



The survival probability is closely related to the MET as we will shortly show. It measures 
the likelihood that, up to time t, the process has been always in the interval [a, £>]: 

S [a , b ] (t - t ; xo) = P{a< X(t) < b, M(t) < b, m(t) > a\X(t ) = x } , 

where we have defined the maximum and the minimum value of X(t), M(t) and m(t), by: 

M(t) = max X(t') , and m(t) = min X(t'). 

t <t'<t t <t'<t 

The financial interest of SP is clear: it may be very useful in risk control. Note, for instance, 
the case b — > oo. The SP measures, not only the probability that you do not loose more 
than a at the end of your investment horizon, like VaR, but also in any previous instant. 

It is notorious that we can recover the MET from the Laplace Transform of the SP. If fact, 
as we have stated above, this is the standard technique used in the literature for obtaining 
METs. The link between both magnitudes becomes apparent if we express the MET in 
terms of P{t[ a ^ < v\xq\, the cumulative distribution function (cdf) of the exit time: 

T[ a ,b]{x ) = / vdP{t [a , b ] < v\x } = / / dudP{t [a>b ] < v\x } 
Jo Jo Jo 

1*00 1*00 l*OQ 

= / dP{t[ ajb ] < v\x }du = / P{t[ a ,b] > u\x }du. 

Jo Jv Jo 

Now, we must realize that the only way that t\ a u can be bigger than any given value is that 
the process has been inside the interval up to that time: 

P {t [a ,b] > t - t \x } = P {a < X(t) < b, M(t) < b, m(t) > a\X(t ) = x } , 

and therefore, 



T[ a ,b](x ) = / S[ atb ](u; x )du = S[ a , b ](s = 0;x ). 
Jo 

It is not surprising that the survival probability follows a renewal equation when also the 
mean exit time can be expressed in such a way — see for instance Ref. [15]. In the present 
case, where we consider that the process properties are depending, at most, on the size of 
last the jump, we can derive the following two-dimensional integral equation for the SP: 

S [aM (t-to;x \AX ) = V(t-t \AX )+ I dt' j dxp{x - x , t' - t \AX )S [aM {t - f; x\AX) 

J to J a 

(7) 



where 

*/£ J — oo 

is the probability that the next sojourn will last more than t — t , given that the previous 
change was of size AX . We can step down the dimension of the integral equation by 
considering the Laplace transform of Eq. (JJJ) , 



S[ ajb ](s; x \AX ) = ^(s|AX ) + / dxp(x - x , s\AX )S [a:b] (s; x\AX). 

J a 

Note that the problem is now much more complex, since it involves the joint pdf of jumps 
and sojourns, p(x,t\AX Q ), not merely its marginal pdfs, h(x\AX ) and ^(r|AXo). Even in 
the fully independent and case, the integral equation is hard to solve: 



S[ a ,b](s;x ) = *(s) + ip(s) / dxh(x - x )S [aM (s;x). 

J a 

The problem of the two-state discrete model without memory, Eq. (J2J), is affordable but 
the complexity of the solution casts few light into the general understanding of the issue. 
Therefore, we have left it for a forthcoming work, and we have focused our attention on the 
symmetric exponential case, Eq. (@J), which gave similar results for the MET — cf Eqs. (jHJ) 
and (jHJ). This model is very suitable for our purposes because reduces the problem from 
solving an integral equation to finding the solution of a second-order (ordinary) differential 
equation: 

dl x S[ a ,b](s; x) = 7 2 (1 -^(s)) S [aM (s; x) - s~ l , 
with the following boundary conditions: 

d x S[ a , b ](s;x = a) = 7 S[ ajb ](s; a) - *(s) , d x S[ a ,b](s; x = b) = -7 S[ ,&](s; b) - *(s) 

Even though, the final expression in the Laplace domain is so intricate: 



S[ a ,b](s; x ) 



l-^(s) 



cosh\jJl -iP(s) (xo- 2 ^) 



1 - ij(s) sinh | 7 ^/l -$(s)L/2 } + cosh <j 7^/1 - i){s)L/2 



that, in general, it cannot be reverted to the time domain. The solution when the process 
begins at the center of the interval is somewhat simpler but still difficult to deal with: 



S[ aib ](s;a+L/2) 



1 - if>(s) sinh | 7X /1 - ijj(s)L/2 } + cosh <j 7^/1 - ip{s)L/2 



The result when the interval width L is infinite, but the process begins at finite distance of 
one of the boundaries, is even shorter: 



S(-oo,x]( s 'i x o) — S{x,oo)(s', x o) — - 

S 



exp < — 7\/ 1 — 4>{s)\x — x \ 



and it can be directly compared with the same outcome for the Wiener process: 



S(-oo,x](s;x ) = S[ X:00 )(s; x ) = - 

S 



1 — exp 



2s, 



\X — Xq\ 



(8) 



where the volatility a is the square root of the diffusion coefficient. The two formulas coincide 
for small values of the Laplace variable s, that is, for large timescales. The resemblance 
between both models when the interval width is bounded is not so evident, because in the 
Wiener case the SP can be only expressed in terms of an expansion series: 



S [aM (s;a + L/2) = 



IL 2 



:-v, 



k=0 



(2k+l)n a 2 n 2 (2k + l) 2 + 2L 2 s' 



In any case, it is easy to check that the long-term behaviour of the MET is similar: 



T [aM {a + L/2)=L 2 /Aa 2 . 



V. RISK CONTROL 

We will finally illustrate how SPs can be used in risk control. In order to clarify the 
concepts we will remove model-dependent inferences by using the outcome corresponding to 
the Wiener case. The Gaussian model is typically used for computing the "Value of Risk" 
(VaR) level. VaR gives the worst return you can obtain at the end of a fixed time interval 
t — to, for a given confidence level a. If we assume that the market volatility is a, then 

VaR = cxv/t^VV- 1 (1 - a) , 

where jV(-) is the cdf for a Normal pdf. This measure of the risk exposure of an open 
position ignores the instantaneous risk aversion of the investor, since it neglects the fact 
that investors may not assume all the paths leading to the same final return. 

This will not be the case if we use SP for quantifying the risk, since it will ensure, within 
the desired level of confidence, that the position is never below the risk measure, which we 
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FIG. 2: Risk values for a position lasting one month, if the process follows a Wiener process with 
a volatility of 15%. We compare, for different confidence levels, the measure of the risk that both 
methods, VaR and SpR, yield. Clearly VaR underestimates the risk. 

will call survival probability risk (SpR). The Laplace inverse transform on Eq. (jHJ) reads, 



In Fig. 121 we will found a comparative example with the two risk measures. 

VI. CONCLUSIONS 

We have argued for the convenience of the use of CTRWs in the modelling of stochastic 
processes in finance. CTRW is a well suited tool for representing market changes at very low 
scales, within the realm of high frequency data. We have shown that this formalism allows 
a thorough description of extreme events under a very general setting: we have obtained 




and therefore 
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renewal integral equations for magnitudes related to these events when the return can be 
described by either an independent or a Makovian process. We have revisited the properties 
of the MET, a statistic that can inform about investment horizons. In previous works we 
found that it seems to scale in a similar way for different assets. We have addressed the topic 
of the SP in finance afterwards. SP has even more severe implications in risk management. 
SpR can improve the efficiency of more traditional methods, like VaR. We have introduced 
new theoretical results on this issue, and shown a practical example of its application. 
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